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Abstract. I give an algorithm for computing the full space of automorphic 
forms for definite unitary groups over Q, and apply this to calculate the au- 
tomorphic forms of level G(Z) and various small weights for an example of a 
^ rank 3 unitary group. This leads to some examples of various types of endo- 
O scopic lifting from automorphic forms for Ui X Ui X Ui and Ui X U2, and to 
1 ^^ ^ an example of a non-endoscopic form of weight (3, 3) corresponding to a fam- 
ily of 3-dimensional irreducible £-adic Galois representations. I also compute 
the 2-adic slopes of some automorphic forms with level structure at 2, giving 
evidence for the local constancy of the slopes. 
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1. Introduction 

If G is a reductive algebraic group over Q, then the spaces of automorphic forms 
for G (of a given level and infinity-type) are well-known to be finite-dimensional 
vector spaces with an action of the Hecke algebra. However, for the vast majority 
of groups it is not known how to actually calculate the dimensions of these spaces 
and how the Hecke operators act. For GL2, the case of classical modular forms, 
there are well-known algorithms based on modular symbols (see e.g. |Ste07aj ) . but 
in general little is known. 

One case in which calculations are possible is where the group G satisfies the 
condition that all arithmetic subgroups are finite. In this case. Gross has shown 
[Gro 99 that the theory of automorphic forms can be set up entirely algebraically, 
without analytic hypotheses. As has been observed by various authors |LP02I 
IDem07j . Gross's spaces are (at least in theory) computable. In this article, I shall 
apply this to an example of a definite unitary group over Q, and show how the 
full space of automorphic forms may be calculated for various small weights and 
levels. We do not need to assume that the class number of G is 1; we do assume 
that the underlying quadratic extension of Q has class number 1, but this is more 
a calculational convenience than a fudamental limitation of the algorithm. 

The results of these computations give rise to some interesting specimens of 
endoscopic liftings of automorphic forms arising from U2 x Ui and Ui x Ui x Ui, 
and also examples of non-endoscopic forms, which give nontrivial families of 3- 
dimensional ^-adic Galois representations. 

The structure of this paper is as follows. Section [2] explains the algorithms 
used to determine the class number of G and the coset decomposition of the Hecke 
operators; these depend on the choice of level group K but not the weight. In 
section |4] we use the results of these computations to calculate the spaces of forms 
of various small weights. In section |5] we shall see how to generalise this to non- 
maximal level groups, and present some examples illustrating the 2-adic continuity 
of the Hecke eigenvalues. In the final section, we discuss some implementation 
details. 

2. Outline of the algorithm 

2.1. Definitions. Recall that if E is an imaginary quadratic number field and 
n > 1, there is a corresponding definite unitary group G — U„,e- which is the 
group scheme over Z defined by 

G{A) ^{ge GL„(A ®z Oe) \ 99^ = 1}, 

where A is an arbitrary commutative ring (and 9^ is the conjugate transpose of g). 
The following properties of G are elementary: 

Proposition 2.1. G(M) is the usual real unitary group [/„, which is connected 
and compact. If p is any prime which splits completely in E, then G is split at p, 
and moreover G(Qp) is isomorphic to GL„(Qp). There are two such isomorphisms, 
corresponding to the two embeddings E ^ Qp, and they are interchanged by the 
inverse transpose involution o/GL„. The centre of G is the norm torus U{1) = 

Since G is compact at infinity, it certainly satisfies the conditions of ^Gro99j : 
so the space of automorphic forms for G of level K (where K is an open compact 
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subgroup of G{Af)) and weight V (where V is an irreducible algebraic represen- 
tation of G over any field F of characteristic 0) is the F-vector space of functions 
/ : G{Af) V such that 

/(75fc) = 70/(5) V7eG(Q), k^K. 

For the rest of this paper, we shall fix n = 3; however, the methods can clearly be 
applied to arbitrary n, and indeed to general totally definite Hermitian spaces over 
Oe (of any integral equivalence class). We shall, however, make the simplifying 
assumption that E has class number 1. Initially, we shall take K = G{Z) = 
Ylg G{Ze), so we obtain the automorphic forms of "level 1" in some sense. 

Remark 2.2. If p is not ramified in E, the subgroup K is a, hyperspecial maximal 
compact subgroup of G{Qp). If p is ramified, then K is special maximal but not 
hyperspecial, unless p = 2 in which case it is not even maximal |Rog90[ §1.10]. 

Theorem 2.3. The irreducible algebraic representations of G over E (or any field 
containing it) are parametrised by triples of integers (a, 6, c) with a,b > 0, where 
the representation Va^.c is the unique highest weight direct summand of Wa,b,c = 
Sym" V (g) Sym'' V* (E) det'^ (where V = Vioo is the standard representation) . The 

central character of Va^b.c is \ z ''+3'=^ 

V y 

Proof. Since G is isomorphic to GL3 over E, this follows from standard results on 
the representation theory of GL3, see e.g. |FH911 §15.5]. □ 

Remark 2.4. We can identify Va^b,c explicitly: there is a natural contraction map 
Wa,b,c Wa-i,b-i,c, and Va,b,c is the kernel of this map. We shall not use this 
remark, however, as in the calculations below we will always have already calculated 
automorphic forms of weight (a — 1,6— 1, c) when we come to deal with those of 
weight (a, b, c); thus we can easily identify forms which we have seen before. 

2.2. Hecke operators and JC-classes. It is well known that the set of double 
cosets 

Gm\GiAf)/K 

is finite. We shall refer to these as iiT-classes. As shown in |Gro99] . if /xi, . . . , /i^ is 
a set of representatives for the isT-classes, then the space of automorphic forms for 
G of level K and weight V is isomorphic to 

r 

where Ti = G{Q) n fj,iKfj,~^, via the map / 1-^ (/(^i), . . . , f{p.r))- Note that the 
groups Tg are arithmetic subgroups of G, and are thus finite. 

If we know a set of class representatives fj,i and the associated groups Ti , then 
it is elementary to calculate the Fi-invariants in Va^b,c for each i and thus read off 
the dimension of the space of automorphic forms. 

The other piece of information we are interested in is the action of the Hecke 
algebra. This is the commutative algebra generated by double cosets KgK, with 
the action of such a coset on automorphic forms being given by 

{[KgK]f) (x)-^/(a;g,) 
3 
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where KgK = \_\gjK. 

So we must decompose KgK into single cosets gjK. For each of these, we could 
choose a coset representative in the form 7j7Zc(j), for some c(j) g {l...r} and 
7j € G'(Q)- This would immediately allow us to calculate (\KgK]f) (1), if / is an 
automorphic form given as an r-tuple /(/^i), . . . fifJ-r), since 



What we actually want, however, is {[KgK]f) (fii); so for each i g {1, . . . , r} we 
need to find elements 7^ G G(Q) and c{i,j) S {1, . . . , r} such that 



i 

For any given element g G G, finding these quantities is a finite search. In fact, 
one has: 

Algorithm A. Let r, s, t he given elements oj G{Kf), which are integral outside 
a finite set S of primes split in E. Then the following algorithm finds all elements 
ofG{Q)r\rKsKt: 

(1) Find some X € Oe such that Xrst, regarded as an element of M^^Af^E) , is 
integral at all places of E. 

(2) Enumerate all matrices g e AI3{Oe) such that gg^ ~ Norm(A). 

(3) For each g in the above list, set 7 = X^^g and calculate the elementary 
factors of r~^jt'^^ at each prime in S. If these coincide with those of s, 
output 7. 

Calculating the elementary factors of a p-adic matrix is easy, but in practice 
one can often short-cut this step, since if VpX is small then there are very few 
possibilities for the elementary factors and these can be distinguished by calculating 
the determinant. The most computationally intensive step is usually (2). For 
convenience, let us define an r-good matrix to be an element of M^^Oe) such that 
gg^ = r. Enumerating all the r-good matrices for a given r is clearly a finite check, 
but possibly a long one if r is large - an algorithm using lattices is described in §6.1| 



As is well known, the Hecke algebra factorises into a restricted tensor product 
of local Hecke algebras at each prime; if p is a prime split in E, and p is a choice 
of prime above p in E, then the local Hecke algebra of at p is generated by the 
double cosets KrjK, where rj is one of the elements 



These are the operators to which we shall apply the above method. We shall 
write Tp i for the operator on automorphic forms corresponding to the double coset 
Krjp iK. Note that Tp_3 is central, and acts on weight {a,b,c) forms as scalar 
multiplication by (p/p)°~''+^^. 

2.3. A refinement. The above discussion assumes that the class representatives 
lii are known. However, we can use a sort of bootstrapping approach to find the 
class representatives at the same time as the Hecke operators. 



{[KgK]f) (1) = ^ /(7,/icO)) = E ° /(^-O))- 



j j 




below. 
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Lemma 2.5. The double cosets Krjp^iK and Krjp 2K both contain + p + I single 
cosets. 

Proof. This is a purely local computation. We shall consider rj = Typ.i! the other 
case is very similar. Let H — G{'Ljy) \ then H acts transitively by right multiplication 
on nonzero row vectors over Fp, and H nrjHrj^^ is the stabilizer of the line (0, 0, *). 
So this subgroup has index |P^(Fp)| = +p+ 1; and thus HrjH should decompose 
into + p + 1 single cosets (which can easily be written down explicitly) . □ 

Now, let us assume that we know a set of representatives for a subset of the class 
set; we can always start with just the principal class G{Q)K, represented by the 
element 1. Let us choose a prime p. Using Algorithm A above, we can find all of 
the single cosets contained in Krjp iK which have a representative of the form 7/i 
with 7 G G{Q) and /i one of the subset of class representatives we know about. 

If we have found less than p'^ +p + 1 cosets, then the lemma shows that we have 
not found the full class set, and moreover, taking a local coset representative gives 
an explicit element fi of G{Af) (supported at p) which is not in any of the _fC-classes 
previously found; so we can add this to the list, calculate the associated group F 
(by applying Algorithm A again, with s = 1) and repeat. 

Clearly, we need a criterion which will allow us to determine when we have found 
the complete class set. To do this, we shall introduce the mass of K, which is a 
quantity strongly related to the class number but in many ways better behaved; 
this can be calculated via the special value of an appropriate L-function. 

Definition 2.6. The mass of a compact open subgroup K of G{Kf), where G is a 
connected reductive algebraic group compact at 00, is the quantity 

Mass(i^) - 

A mass formula for unitary groups is given in |GHY01| . This applies to a some- 
what more general class of unitary groups than we consider here, so we shall state 
a special case: if G is the unitary group C/„,_e as defined above and n is odd, and 
K = G{Z), then 

Mass(if) = ^.L(M).r(G)-nAp. 

p 

Here L{M) is the value at s = of the L-function of the motive of G; t(G) is the 
Tamagawa number of G; and the local factors Ap are 1 unless p is ramified in _B, in 
which case they are i. (For n even, there is a similar formula but the local factors 
are all 1.) 

Formulae for t(G) and L{M) are known, and are given in |GHYOl] . For unitary 
groups we have t(G) = 2; and if x is the Dirichlet character corresponding to the 
quadratic extension i?/Q, then 

n 

L{M) = \{L{\-t,x1. 

r=l 

which is easy to compute using generalised Bernoulli numbers as in chapter 4 of 
|Was82j . 

Remark 2.7. Knowing the mass in advance allows one to determine whether the 
full class set has been found; but it is not completely clear that this will ever 
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actually occur - that is, that every K -class occurs in the decomposition of some 
Hecke operator. On the other hand, one can show (using the assumption that E 
has unique factorisation and strong approximation for the derived subgroup S'J/s.e, 
cf [PR94, ch. 7]) that every ii'-class does have a representative supported at p for 
every prime p. Hence if computing Tp.i and Tp.2 for enough p does not find enough 
if-classes, one can switch to computing multiple Hecke operators at some fixed 
small prime. However, I am not aware of any case in which this is necessary. 

3. Example: E ^ Q{y^) 

To illustrate the algorithm above, we shall carry out the calculations in the case 
E — Q{V—7). This choice is convenient, as then G is split at 2, and — ±1, so 
G(Q) K is small (it is the group of monomial matrices with entries in O^, and 
thus has order 48). Let w = , so Oe = Z[a;]. 

We shall begin by calculating the mass of K. Here x is the Kronecker symbol 
(=^); using standard methods we find that L{0,x) = 1, L{—l,x^) = — and 
i(-2,x^) = so L{M) = Hence 

14 11 

Mass(is:) = 2 - - = — . 

^ ^ 8 21 2 42 

Since > j^, the mass of the principal X-class, there must be other iiT-classes; 
we shall find them by the bootstrapping argument explained above. 

If p is any split prime and p a choice of factor of p in E, then we can apply 
Algorithm A to find the set G{Q) nKij^^iK, taking A = p. For instance, let us take 
p = 2 and choose p to be the ideal generated by oj. Then we find 288 such matrices; 
since G(Q) n K has order 48, this gives us 6 right cosets 'jK, which are represented 
by the three diagonal matrices with entries 1, 1, in some order, and the three 
permutations of the rows of 

^10 

2 
vO 2 

Lemma |2.5| tells us that there is one coset missing, and the proof of the lemma 
constructs an element of this coset, namely the element a of G(A/) that is 1 at all 
primes p 2 and under the isomorphism G(Q2) — GL3(Q2) corresponds to 

1 0^ 
1 

-1 -1 2^ 

The associated group Fq, has order 336; indeed it is isomorphic to the direct 
product of ±1 with a simple subgroup isomorphic to PSL2{V^). Conveniently, this 
group contains G(Z) as a subgroup of index 7; it is generated by G(Z) and the 
element 

1 l-UJ 

' 2 2 

1 --l+oj 

2 2 

- 
N 2 " 

of order 7. 

Since ^ + = = Mass(_ft'), that's all: the class number of G is 2. 

We can now complete the calculation of the Hecke operator Tp.i as follows. In 
the notation of the previous section, we are taking fj,i = 1 and fi2 — ol. We have 
decomposed ii\Kr\K as a union of 6 cosets of the form jfiiK and 1 coset 7/i2-R^. 
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Searching for rational elements of ii2KriK (using Algorithm A with A = (1 — loY) 
finds 7 cosets already, so we do not need to consider the fourth possibility, which 
is whether there are any elements of G'(Q) n fi2Kr]Kii2^ ■ 

Since we have now found the complete class set, it is clear how to calculate Tp ^ 
for any split prime p. Note that we need to apply Algorithm A four times, with 
A = p, ujp, ujp and 2p, so we need to enumerate r-good matrices for r — p, 2p and 
4p. 

Remark 3.1. Calculating Hecke operators at the inert primes is also possible, but 
more awkward both theoretically and computationally. Since the maximal split 
torus in G{Qp) for p inert has rank 1, there is only one Hecke operator to deal with, 
but calculating its coset decomposition requires enumerating r-good matrices for 
r — p^ , which rapidly becomes computationally infeasible. Also, the eigenvalues of 
these Hecke operators are rather more indirectly related to Galois representations. 

4. Experimental results: level 1 

In this chapter I'll summarise the calculations I have made for various small 
weights. Note that there are no forms of weight (a, 6, c) if a + & + c is odd since the 
level group contains —1; and changing c by a multiple of 2 just induces a twist, so 
I shall restrict to the case where c = if a + 6 is even and c = 1 if o + & is odd; this 
is what is meant by "weight (a, 6)". 

4.1. Dimensions of spaces. These can be deduced from classical results of in- 
variant theory for actions of groups on polynomial algebras. See table [l] 







a 









1 


2 


3 


4 


5 







2 





1 





3 







1 











1 


1 


2 


b 


2 


1 





2 


1 


5 


2 




3 





1 


1 


4 


3 


6 




4 


3 


1 


5 


3 


8 


7 




5 





2 


2 


6 


7 


10 



Table 1 . Dimension of spaces of automorphic forms 



Note that if a, h are small and a + 6 is odd, then the spaces are often zero. Thus 
I will mostly restrict to the case of even weight. 

4.2. Forms arising from Grossencharacters. For small weights, many of the 
forms that arise appear to have Tp_i eigenvalues of the form Xi(p) + X2(p) + 
X3(p), where Xi{p) — P'^P" for some (r, s). Note that there exists an unramified 
Grossencharacter of E sending a uniformiser at p to p'^p'' whenever r-|- s is even, as 
E has unique factorisation and unit group ±1; let us write x(r, s) for this character. 

In all cases other than the constant function, we find that the characters arising 
have r + .s = 2. In table [2] below I list these, based on calculating Hecke operators 
at all split primes p < 60. (I have assigned arbitrary numbers to Gal(£^/i?)-orbits 
of eigenforms, in order of discovery.) 
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Weight 


Form 


Character 


(0,0) 


1 


x(o,o)- 


f x(i,i) + 


X(2,2) 


(0,0) 


2 


X(2,0)- 


f x(i,i) + 


X(0,2) 


(2,0) 


1 


x(4,-2) 


+ x(i,i)- 


HX(0,2) 


(4,0) 


3 


x(6,-4) 


+ x(i,i)- 


I-X(0,2) 


(2,2) 


2 


x(4,-2)^ 


f x(i,i) + 


x(-2,4) 



Table 2. Forms arising from Grossencharacters 



4.3. Forms arising from classical modular forms. Other forms of small weight 
appeared to have Tp^i eigenvalues given by x(^i, ^i) + x(''2, S2)ap for ap the Hecke 
eigenvalue attached to a classical modular newform. Note that this class includes 
all of the forms above (some in multiple ways), since for each r there exists a CM- 
type modular form of weight r — 1 with Op = x(r, — r) + x(— r). The automorphic 
forms on U{3) that are associated with non-CM-type forms are listed in Table [3] 
below. 



Weight 


Form 


Character 


Modular weight 


Modular label 


(4,0) 


1 


X(l,l)+x(0,-4)ap 


7 


7k7B[3]l 


(4,0) 


2 


x(0,2)+x(l,-4)6p 


6 


49k6El 


(3,1) 


1 


X(l,l)+x(-l,-3)ap 


7 


7k7B[3]l 


(2,2) 


1 


x(l,l)+x(-2,-2)ap 


7 


7k7B[3]l 


(3,3) 


1 


x(l,l)+x(-3,-3)cp 


9 


7k9B[3]l 



Table 3. Forms arising from non-CM-type modular newforms 



In the third column of the table, I have reused letters when the same modular 
form occurs more than once. The labels in the final column are those used in 
William Stein's online tables. 

Note that above x(r, s) was defined only for r + s even. If r + s is odd, there is 
no unramified Grossencharacter sending a uniformiser at p to p'"p'^ - indeed, this 
is not well-defined as the choice of generator of the ideal p is only defined up to 
±1 - but there is a unique Grossencharacter which is unramified outside A = V—7, 
whose restriction to 0{E\) is the order 2 character of , and whose restriction to 
Eao = C is z 1-^ z~^z~'^; this character sends a uniformiser at p to (^) p^p'', which 
is independent of the choice of generator of p. This is what is meant by x(l, — 4) in 
the above table. 

4.4. Interpretation: Satake parameters. Recall that the Satake parameters of 
an unramified representation tt^ of GL„(Qp) are the eigenvalues of the semisimple 
conjugacy class that is the Langlands parameter of the representation TTp. These 
can be obtained as the roots of the Satake polynomial 

n 
i=0 

where tp i is the eigenvalue of Tp i on the GL3(Zp)-invariants of Tip] see |Che03[ 
§1.8.4]. 

We introduce the following definition, from |BC07i §6.4.5]: 
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Definition 4.1. Let X be the multiset of Satake parameters (with multiplicities); 
then TTp is is said to be almost tempered if there is a partition of X into subsets Xi 
such that each Xi is of the form {y, y/q, . . . , y/?'"'^^}, for rrii = \Xi\, and the real 
number \ Yixex a^l^^™' does not depend on i; here \ ■ \ denotes the complex absolute 
value. 

If TTp is tempered, which means in ths context that it is a subquotient of the 
parabohc induction of a unitary character of the standard Borel subgroup, then all 
of the Satake parameters have complex absolute value and so it is almost 

tempered with all the Xi singletons; but not all almost tempered representations are 
of this type - for example, a representation with Satake parameters {p^^^ ,p~^^'^ , 1} 
would be almost tempered^ However, if TTp is generic (admits a Whittaker model) 
and almost tempered, it must be tempered. 

A conjecture due to Arthur predicts that if TTp is the local factor of a cuspidal 
automorphic representation of any unitary group, then TTp is almost tempered, which 
implies that if Wp is generic, all the Satake parameters must have complex absolute 
value p. See |BC07[ §6.4] for further discussion. 

Since our form of U^^e is compact at oo, all automorphic forms are cuspidal. For 
all the forms we have encountered, the local factors at all split primes calculated 
are generic (except for the constant form, which is never generic). For the forms in 
Table [2] and Table [3] above, the Satake parameters are easy to read off: for Table 
|2] the Satake parameters are the values at p of the three characters listed, and for 
Table [3] the Satake parameters for xi + X20.p are {xi(p)i X2(p)a, X2(p)/3} where a 
and P are the Satake parameters associated to the form /. The normalisation that 
is conventional for classical modular forms of weight k gives both a and f3 absolute 
value pC^"!)/^. Thus the characters x(r, s) in Table [2] are forced to have r + s = 2; 
and for characters of the form x(''ii ■Si) + x(''2, S2)ap, we'd better have ri + si — 2 
and r2 + S2 = 3 — k. 

4.5. Galois representations. Let / be an eigenform of weight (a, b, c) and level 
G(Z). By results of Blasius and Rogawski [BR92j . for each prime £ there is a 
semi-simple representation 

Pf : Gal(E/E) GL3(Q,) 

unramified outside 7£ and crystalline at i, such that for p — pp any other prime 
split in E not dividing 7i the characteristic polynomial of p/(Frobp) should be the 
Satake polynomial of / at p as defined above. 

(The proofs of this statement in the literature assume that the weight is "reg- 
ular", i.e. a ^ and b ^ 0. For non-regular weights, ^-adic interpolation in the 
manner of |Che03| Ch. 1] gives a representation with values in GL3(C£) as long as 
£ splits in the field E, but it is not clear whether it is crystalline. It is, however, 
necessarily Hodge-Tate.) 

Here tp^i is the eigenvalue of Tp^i on /, so in particular tp^i is the trace of Frobe- 
nius. Thus if the formulae conjectured above for the Tp^i eigenvalues of small weight 
forms are correct, the associated representations pf decompose as direct sums of 



In general, the almost tempered representations are exactly those that are twists by unramified 
characters of parabolic inductions from a unitary character of the Levi subgroup of some parabohc; 
if this parabolic is not minimal, the representation will not be tempered. The above example 
corresponds to the trivial character of the Levi subgroup GL2 X GLi C GL3. 
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smaller-dimensional representations (sums of three characters, or of a character and 
a twist of a 2-dimensional modular representation) . 

If f is a prime split in E, then we can restrict the representation to a decom- 
position group at either of the primes above £. Since this local representation 
is crystalline, it is certainly Hodge- Tate; but the weights at these two different 
primes are not generally the same - for one decomposition group the weights will 
be —b+c, c+ 1, a + c+2 and for the other the weights will be —a — c, —c+1, b — c+2. 
This fits the tables above, as x(a, b) has weight a at one prime and b at the other, 
and the representation associated to a modular form of weight k has weights and 
k — 1 (whichever prime is used) . 



4.6. Matching weights. Given a fixed weight (a, 6, c) G N x N x Z, the constraints 
of the valuations of Satake parameters and the Hodge- Tate weights at the two 
primes above £ leave few possibilities for endoscopic forms at this weight. 

Firstly, let us consider the case of automorphic forms / arising from the endo- 
scopic subgroup Ui X Ui X Ui. These correspond to Hecke eigenvalues that are the 
sum of three Grossencharacters of E, as in table |2] If / is to have level 1, then the 
three Grossencharacters must be unramified; so the Tp^i eigenvalues must be of the 
form x(''i) si) + x(''2, S2) + x(''3, S3) with -|- Sj even for each i. 

We shall suppose c — 0. Then if the local factor of / is tempered at some prime 
split in E, we must have ri + Si — 2 for each i. Then it is easy to see that the only 
possibility is x(a -|- 2, —a) -I- x(l, 1) + x{~b, b + 2), which has level 1 if a and b are 
both even. This explains why we have seen one eigenform of this type when a, b 
are both even and none otherwise. 

Assuming Arthur's conjecture, if / is not tempered, then it is almost tempered, 
two of its Satake parameters have ratio p, and we must have (up to a permutation) 
either ri + si = 4, r2 + S2 = 2, and ^3 -I- S3 = 0, or ri + si = 2, r2 + S2 — 3, and 
''a + S3 = 1. In the first case we also must have ?'2 = ^1 ^ 1 ^md r3 = ri — 2, and 
the same for the s^; this can only occur for the constant form of weight (0, 0). The 
second case cannot occur at level 1, as we must have -I- Si even. 

The situation for endoscopic lifts from U2 x Ui is more complicated. Let us 
assume our putative automorphic form / on U3 is generic at all places, has weight 
(a, 6,0) for integers a,b > 0, and has Hecke eigenvalues of the form xi^i^^i) + 
x(^2, S2)ap where Op is the Hecke eigenvalue of a classical modular cusp form g of 
weight k. (Note that this includes the previous case, as endoscopic lifts from Ui x Ui 
to U2 correspond to modular forms g of CM type.) The Galois representation 
associated to g has Hodge-Tate weights and k — 1, and its Satake parameters 
at p have complex absolute values pC^"!)/^, so we must have ri -I- si —2 and 
''2 + S2 = 3 — fc. Checking all the possible cases, we find that there are potentially 3 
families of automorphic eigenforms on U3, with Tp 1 eigenvalues and weights given 
by 

x(l, 1) + x(^s, 3 — fc -I- s)ap weight (fc — 3 — s, s, 0) 
X(2 -r,r) + x{l,2- k)ap weight (fc - 2, r - 2, 0) 
X(r, 2 - r) -I- x(2 - fc, l)ap weight (r - 2, fc - 2, 0). 
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4.7. A non-endoscopic form. Weight (3, 3) is the first case where the action of 
the Hecke operators is not diagonalisable over E. One obtains two Hecke-invariant 
subspaces of dimension 2, corresponding to Galois orbits of eigenforms. 

The first subspace sphts over £'(\/46), and the Tp.i -eigenvalues are 1) + 
x(— 3, — 3)cp, where Cp is the unique Galois orbit of modular eigenforms of weight 
9 for ri(7) such that C2 is defined over this field. In this case, Tp 2 coincides with 

The second space is more interesting. Here Tp 1 and Tp 2 commute with each 
other and have the same eigenvalues, but are not identical; so for the two eigenforms 
in the orbit, the Tp i-eigenvalue of one form is the Tp,2-eigenvalue of the other form 
and vice versa. The eigenvalues are all defined over Q(-\/— 259), and the Satake 
polynomial for one of the primes above 2 is 

Since the Satake polynomial is irreducible over 259), the Satake parame- 

ters are all Gal(i?/i?)-conjugate; so none of them arise from Grossencharacters of 
E, and thus the corresponding ^-adic Galois representations must be irreducible of 
dimension 3. 

Let's check the valuations of the Satake parameters. It can be shown |Loe08] 
that a cubic over C has all its roots of equal valuation if and only if, after scaling 
so the constant term is 1, it is of the form + aT^ + aT + 1 with a lying inside a 
deltoid in the complex plane with centre and one cusp at 3; this is equivalent to 
a — X + iy with x, y real such that 

C{x, y) = (a;2 + y'^f - Sx{x^ - 3y^) + I8{x^ + y^) - 27 < 0. 

In this case we find that C{x, y) — — 4^9/ , so the condition is satisfied. 



5. Some higher level examples 

We shall now introduce some level structure at the prime 2. In preparation for 
future work in which I intend to consider overconvergent automorphic forms, we 
shall use a slightly different definition of automorphic forms in which the action 
is twisted to the other side. We shall fix a prime p which is split in E and an 
irreducible algebraic representation of G (which will be defined over Qp), and 
consider functions 

f:G{Kf)^V such that .f{l9k) = k-\f{g) 

for all 7 S G(Q) and k & K ion some open compact K. This is isomorphic to the 
space considered above, via the map f{g) ^ gp°f{g) (if we fix a choice of embedding 
oi E Qp), but has the advantage that we only ever need to consider the action 
on V of elements of Kp] this is useful when considering p-adic interpolation, as in 
|Che03j or my University of London PhD thesis (in preparation) . 

The action of the Hecke algebra on this space is slightly different: we define 
an action of a double coset KgK on automorphic forms of level K by writing 
KgK = \_\giK, and defining 

{[KgK]o f){x) ^Y.^g,)pf{xg,). 
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It's easily verified that this gives a well-defined linear endomorphism of the space 
of automorphic forms. 



5.1. X-classes. We shall set p = 2; and wc will define our level group L to be 
G{Ze) at all ij^ 2, and at 2 the subgroup Lp of GL3(Z2) which reduces mod 2 to 

(* * *\ 
0**1. 
* */ 

This is the parahoric subgroup associated to a parabolic subgroup of GL3. We now 
need to calculate L-classes and decompositions of Hecke operators for this case. We 
shall do this by comparing L to the group K = GL3(Z) considered above, using 
the following lemma: 

Lemma 5.1. Let U' C U be any finite index subgroups of G{Af). Then we have 
Mass(C/') = [U : U'] Mass{U). 

Proof. This is an essentially trivial but fiddly manipulation. 



Mass(?7') = X ^ 



\G(Q)nvU'v-^ 



E E 



\G(Q) n ukU'k-^n-^ size of orbit of k 

(where we consider the action of f/ fl /Lt~^G(Q)/ti on U/U' by left multiphcation) 

_ ^ ^ 1 \G{Q)ni^kU'k-^ii-^\ 

~ 2^ 2^ \G(q)niikU'k-^ij-M \G(Q)nijUi2-^ 

[lJ.]eG{Q)\G{Af)/U keU/U' I ^ ^1 \ \-^J M' r- \ 

= Mass(C/) • I E ^ 
\keu/u' 

= Mass(f/) • [U : U'] 

□ 



It's readily seen that Lp is the subgroup of matrices in GL3(Z2) whose reduction 
mod 2 stabilises (*,0,0)'^ under the left multiplication action on column vectors; 
GL3(Z2) clearly acts transitively on these, so [K : L] = |P^(F2)| = 7. Hence 
Mass(L) = ^ = |. 

The proof of the above lemma indicates how to work out the L-classes, by cal- 
culating, for each ii'-class G{Q)fiK, how K C] n~^G{Q)fi acts by left multiplication 
on K/L, which we can identify with P^(F2). 

For the identity ii'-class, the group G(Q) fl K consists of a product of diagonal 
±l's (which all reduce to the identity and can be ignored) and permutation matrices. 
So the action corresponds to permutation of coordinates on P^(F2), which has 
exactly 3 orbits classified by the number of nonzero terms. Coset representatives 
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are given by: the identity (having stabiHser of order 16); the matrix 110 



For the nonidentity _ft'-class, the extra generator r of G{Q) n aKa ^ conjugates 
to an element of G{Af) whose matrix reduces mod p to 



This acts transitively on the elements of ¥'^{¥2), so G{Q)aK = G{Q)aL. The 
stabilizer group must thus have size 48, and it is easy to enumerate its elements as 
it is a subgroup of the finite group Fq, of the previous two chapters. 

5.2. The operator U. Motivated again by the theory of overconvergent auto- 
morphic forms, we consider the action of the Hecke operator corresponding to the 
double coset Lirjp 2L, where p is the prime dividing 2 generated by oj. As the role it 
plays is roughly analogous to the Atkin-Fehner J/p-operator in Coleman's theory, 
we shall denote it by the letter U . 

To decompose L?7p_2-^ into single cosets locally, we observe that if 77 = '>lp,2, 
then Lp n r]Lpr]~^ consists of those elements of Lp whose first column reduces to 
(*,0,0)-^ mod 4. This is an index 4 subgroup of Lp and coset representatives 
may be obtained by taking the first column to be (1,0, 0)^, (1, 2, 0)"^, (1, 0, 2)^ and 
(1, 2, 2)-^. We hence obtain a decomposition of LrjL as |Ji=i 't'iVL- 

As in the level K case, the final step is to write each single coset in terms of our 
chosen L-class representatives ^1 ... ^4, and more generally to do the same for the 
products firff'sV^- This is easy using an adaptation of Algorithm A: we note that 
if 7 G G{Q) satisfies 7/it S iir4'sVL, then the largest denominator that can possibly 
occur in 7 has norm either 2, 4 or 8, depending on (r, s, t). We then generate lists 
of the 2, 4 and 8-good matrices, and for each coset fir^'s'^L, we use these to find 
a representative for this coset in the form 7/^,4 (testing each of the possible values 
of t in turn until one works). Note that we must explicitly test (pj^ fi~^jfj,^ for 
membership of L, rather than dodging this check using elementary factors as before. 
The results of this computation are shown in Table |4j 

Remark 5.2. Since the group Lp we have chosen is a parahoric subgroup, one can 
attach invariants to double cosets LpgLp C GF3(Qp) (analogous to elementary 
factors) using the Bruhat-Iwahori decomposition. One can probably use this to 
calculate the Hecke operators more efficiently; but the slower method above is 
more flexible and would apply to any group L for which we could compute K/L 
and L/{L n riLrj^^). 

5.3. Automorphic forms and slopes. The coset representatives calculated in 
Table |4] above can now be used to calculate automorphic forms of level L using the 
same machinery as before. The choice of level structure and of U operator should 
imply that the resulting objects interpolate well 2-adically as the a component of 
the weight varies. See Table [5] for a table of the 10 smallest 2-adic slopes of the U 



1 



(again with 16); and the matrix 1 1 with stabiliser of order 48. 
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1^3 L 



Table 4. Decomposition of the double coset LrjL 



a 


Slopes 



2 
4 
6 
8 
10 
12 
14 
16 


12 2 

1 1 1 3 4 

012||334 5 5 

ni33qq7 7 11 11 

"-^22''''2 2 2 2 

0122244| 1 1 
0l||3555 5 f 
012|||34 5 6 

r)i33oo7 7 11 11 
^^22 2 2 2 2 
0122244^f^ 



Table 5. The 10 smallest slopes of U acting on forms of level L 
and weight (a, 0,0), for even a < 16. 



operator on forms of weight (o, 0, 0), for various integers a (necessarily even). This 
table shows clear evidence of the 2-adic local constancy of the slopes. 

Note that the slopes of U are not remotely locally constant with regard to vari- 
ation in the b direction in the weight lattice; indeed the dimension of its ordinary 
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subspace grows as b gets large (in the archimedean sense!). This is unsurprising, as 
U does not define a compact endomorphism on the space of overconvergent forms 
as defined in j Che03| . but it is compact restricted to the smaller space of "semi- 
classical" automorphic forms associated to the parabolic subgroup corresponding 
to the parahoric Lp - see my forthcoming University of London PhD thesis for 
the definitions and properties of these spaces. This space corresponds to allowing 
p-adic variation in the a direction but not the b direction in the weight lattice. 

A striking aspect of these tables of slopes is the fact that they are often non- 
integral. This is very much contrary to the case of GL2 where slopes are "usually" 
integral, at least for weights in the centre of weight space (Wan has conjectured 
that for any fixed weight k the denominators of the slopes of overconvergent forms 
at weight k is bounded) . Extending the computations above gave an example of a 
form with slope at weight (5, 5, 0). 

6. Implementation details 

6.1. Algorithms for r-good matrices. The time-consuming step in Algorithm 
A, which was used in several places, is to calculate all elements of M3{Oe) satisfying 
mm^ — r; to calculate the Hecke operator at a split prime p one needs to consider 
r — p, 2p, and 4p, so to have a reasonable supply of primes one needs to consider 
r « 100. Clearly this is equivalent to finding all orthogonal triples of vectors 
of length r in 0|,, regarded as a 6-dimensional Z-lattice. I experimented with two 
different algorithms, both of which have the same first step of generating a complete 
list £ of the vectors of length r. This list is typically quite large - since the theta- 
function of the lattice is a modular form of weight 3 and level ri(7), it is easily 
seen that for r = p, 2p,Ap with p a prime split in E, there are O(p^) such vectors 
(with leading terms that can be calculated precisely using Eisenstein series in each 
case). 

The obvious algorithm is to regard this as a graph triangle enumeration prob- 
lem: consider the graph whose vertices are the vectors in C and where vertices 
u, V are joined by an edge if they are orthogonal. This graph is typically not very 
dense, so from an appropriate sparse representation of its adjacency matrix one can 
enumerate all triangles in it quickly by a vertex-iterator approach. The most com- 
putationally difficult step is calculating the adjacency matrix: testing all possible 
pairs of elements of £ will take 0{p*) steps. 

A more devious alternative is to calculate, for each vector v of length r, a basis 
for the orthogonal complement of v in O^, as a 4-dimensional Z-lattice; enumerate 
vectors of length r in this; and find the orthogonal complement of these, to give 
orthogonal triples. This option turns out to be very much quicker for large r, as 
the orthogonal complement lattices typically have very few short vectors. 

6.2. Programs included with this paper. A selection of computer programs 
that were used in the calculations above can be downloaded from my personal 
website at http : //www. ma. ic . ac .uk/~dl505/maths/progrcmis/programs .html 

These use a combination of the Sage |Ste07bj and Magma p3CP97! computer alge- 
bra systems; users who do not have access to Magma can still use the code for level 
1 automorphic forms, but will be restricted to those primes for which precalculated 
lists of r-good matrices are included. 



16 



DAVID LOEFFLER 



7. Acknowledgements 

I would like to thank Gaetan Chenevier, for initially raising the question of 
computing automorphic forms on unitary groups and for many helpful discussions 
since, and my PhD supervisor Kevin Buzzard, for his unfailing encouragement and 
support. 

I am also grateful to Michael StoU, for suggesting to me the second, more efficient 
lattice-based algorithm for r-good matrices described in §6. 1| above. 

References 

[BC07] Joel Bellaiche and Gaetan Chenevier, p-adic families of Galois representations and 

higher rank Selmer groups, Asterisque (to appear), 2007. 
[BR92] Don Blasius and Jonathan D. Rogawski, Tate classes and arithmetic quotients of the 

two-ball, The zeta functions of Picard modular surfaces, Univ. Montreal, 1992, pp. 421- 

444. M R 1155236.1 

[BCP97] Wieb Bosma, John Cannon, and Catherine Playoust, The Magma algebra system, I. 
The user language, J. Symbolic Comput. 24 (1997), no. 3-4, 235-265. MR 1484478. 

[Chc03] Gaetan Chenevier, Families p-adiques de formes automorphes et applications aux con- 
jectures de Bloch-Kato, Ph.D. thesis, Universite Paris VII, 6 2003. 

[Dem07] Lassina Dembele, Quaternionic Manin symbols, Brandt matrices, and Hilbert modular 
forms. Math. Comp. 76 (2007), no. 258, 1039-1057 (electronic). MR 2291849. 

[FH91] William Fulton and Joe Harris, Representation theory: a first course. Graduate Texts 
in Mathematics, vol. 129, Springer, 1991. MR 1153249. 

[GHYOl] Wee Teck Gan, Jonathan P. Hanke, and Jiu-Kang Yu, On an exact mass formula of 
Shimura, Duke Math. J. 107 (2001), no. 1, 103-133. ,MR 1815252., 

[Gro99] Benedict H. Gross, Algebraic modular forms, Israel J. Math. 113 (1999), 61-93. |MR| 
1729443. 

[LP02] Joshua Lansky and David Pollack, Hecke algebras and automorphic forms. Compos. 

Math. 130 (2002), no. 1, 21-48. MR 18836907 
[Loe08] David LoefHer, Adventures with polynomials: a criterion for Weil numbers. Eureka 

59 (to appear) (2008), Available from ht tp: //www.ma. ic ■ ac.uk/~dl505/maths/cublcs ■ | 

pf 

[PR94] Vladimir Platonov and Andrei Rapinchuk, Algebraic groups and number theory. Pure 
and Applied Mathematics, vol. 139, Academic Press, 1994. MR 1278263. 

[Rog90] Jonathan D. Rogawski, Automorphic representations of unitary groups in three vari- 
ables, Annals of Mathematics Studies, vol. 123, Princeton Univ. Press, 1990. MRj 
[r081540.l 

[Ste07a] William Stein, Modular forms, a computational approach. Graduate Studies in Mathe- 
matics, vol. 79, Amer. Math. Soc, 2007. MR 2289048. 

[Ste07b] , Sage Mathematics Software (Version 2.8.9), The Sage Group, 2007, [http : | 

//www . sagemath . org/ 

[Was82] Lawrence C. Washington, Introduction to cyclotomic fields, Graduate Texts in Mathe- 
matics, vol. 83, Springer, 1982. M R 1421575.. 
E-mail address: david.loeffler050imperial.ac.uk 

Department of Mathematics, Imperial College, South Kensington, London SW7 2AZ 



